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Relations between critical exponents, or scaling laws, at both continuous and discontinuous quan¬ 
tum phase transitions are derived and discussed. In general there are multiple dynamical exponents 
at these transitions, which complicates the scaling description. Some rigorous inequalities are de¬ 
rived, and the conditions needed for these inequalities to be equalities are discussed. Scaling laws 
involving the specific-heat exponents that are specific to quantum phase transitions are derived and 
and contrasted with their counterparts at classical phase transitions. We also generalize the ideas 
of Fisher and Berker and others for applying (finite-size) scaling theory near a classical first-order 
transition to the quantum case. We then apply and illustrate all of these ideas by using the quan¬ 
tum ferromagnetic phase transition in metals as an explicit example. This transition is known to 
have multiple dynamical scaling exponents, and in general it is discontinuous in clean systems, but 
continuous in disordered ones. Furthermore, it displays many experimentally relevant crossover 
phenomena that can be described in terms of fixed points, originally discussed by Hertz, that ulti¬ 
mately become unstable asymptotically close to the transition and give way to the asymptotic fixed 
points. These fixed points provide a rich environment for illustrating the general scaling concepts 
and exponent relations. We also discuss the quantum-wing critical point at the tips of the tricritical 
wings associated with the discontinuous quantum ferromagnetic transition from a scaling point of 
view. 

PACS numbers: 05.30.Rt; 05.70.Jk; 75.40.-s 


I. INTRODUCTION 


Classical or thermal phase transitions are special 
points in parameter space where the free energy is not 
analytic as a consequence of strong thermal fluctua¬ 
tions, see, e.g., Refs. mm It is customary to distin¬ 
guish between second-order or continuous phase tran¬ 
sitions, where the order parameter (OPP goes to zero 
continuously at the transition, and first-order or discon¬ 
tinuous ones, where the order parameter displays a dis¬ 
continuity. Classic examples for the former are the fer¬ 
romagnetic transition at the Curie temperature in zero 
magnetic field, or the liquid-gas transition at the criti¬ 
cal point; for the latter, a ferromagnet below the Curie 
temperature in a magnetic field, or the liquid-gas tran¬ 
sition below the critical pressure. The nonanalytic free 
energy translates into singular behavior of observables. 
At a second-order transition, this usually takes the form 
of power laws, in which case the singularity is character¬ 
ized by critical exponents. Let m be the order param¬ 
eter (in the case of a ferromagnet, to is the magnetiza¬ 
tion, in the case of a liquid-gas transition, the density 
difference between the phases), h the external field con¬ 
jugate to the order parameter, and r the dimensionless 
distance from the critical point for h = 0 ® at a thermal 
phase transition, r is usually chosen to be the distance 


from the critical temperature T c , r = T/T c — 1. Then 
the specific heat C behaves as C oc |r|~ a , the order pa¬ 
rameter to vanishes according to m(r —4 0—) oc (— r)@, 
the order-parameter susceptibility \ diverges according 
to x{T) oc |r| -7 , etcP^a, /3, and 7 are examples of critical 
exponents. Underlying all of these singularities is a di¬ 
verging length scale, the correlation length £. It measures 
the distance over which order-parameter fluctuations are 
correlated, and diverges as f oc |r| _ ", which defines the 
exponent v. Two other important critical exponents are 
6 , which describes the behavior of the order parameter 
as a function of its conjugate field h at criticality, m(r = 
0 , h) oc /i 1 /" 5 , and the exponent 77 , which describes the 
decay of the order-parameter correlation function at crit¬ 
icality, G(r = 0, |ao| —¥ 00 ) = (m(x)m( 0)) oc l/\x\ d ~ 2+ri , 
where d is the spatial dimensionality. 


The critical exponents are not all independent. At 
most classical transitions, it turns out that specifying two 
exponents determines all of the others. This was recog¬ 
nized early on and put in the form of exponent relations 
(also referred to as “scaling relations”, or “scaling laws”; 
the latter not to be confused with the homogeneity laws 
that are often called “scaling laws”). For instance, the six 
exponents a, /3, 7 , 5 , rj , and v are related by the four 
relations listed in Eqs. (B1 1 . These were initially derived 
from scaling assumptions, i.e., generalized homogeneity 
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laws for various observables near a critical point, and 
later understood more deeply in terms of the renormal¬ 
ization group (RG), which allowed for a derivation of the 
homogeneity laws. Some of the exponent relations are 
more robust than others. For instance, for a ^> 4 -theory 
in d, > 4 spatial dimensions, the Widom, Fisher, and 
Essam-Fisher scaling relations expressed in Eqs. 
main valid, whereas the hyperscaling relation, Eq ), 
fails. This is related to the notion of dangerous irrele¬ 
vant variables (DIVs), i.e., coupling constants that flow 
to zero under RG transformations but that some observ¬ 
ables depend on in a singular way. As a result, “strong 
scaling”, i.e., a simple homogeneity law for the free en¬ 
ergy, breaks down. More generally, there are constraints 
on the critical exponents that take the form of inequali¬ 
ties. They rely on much weaker assumptions than strong 
scaling and in some cases are rigorous. Examples are the 
Rushbrooke inequality, Eq. ( ]B5| ) , and the lower bound for 
the correlation-length exponent v in disordered systems 
that is discussed in Appendix |C| 

A quantum phase transition (QPT) occurs by defini¬ 
tion at zero temperature, T = 0, as a function of some 
non-thermal control parameter such as pressure, compo¬ 
sition, or an external magnetic field.® The critical behav¬ 
ior at T = 0 is governed by quantum fluctuations rather 
than thermal ones. The role of temperature is thus dif¬ 
ferent from that at a thermal transition, which neces¬ 
sitates the introduction of additional critical exponents. 
For instance, the order-parameter susceptibility will vary 
as |r| -7 , with r the dimensionless distance at T = h = 0 
from the critical point, but as T _7T as a function of tem¬ 
perature for r —> 0, with xt in general different from 
7 . Similarly, we need to define exponents and ut 
that describe the behavior of the order-parameter and 
the correlation length, respectively, as functions of the 
temperature in addition to r. 

The case of the specific heat c is less straightforward 
since c vanishes at T = 0 even away from any critical 
point. Related to this, in the thermodynamic identity 
that underlies the Rushbrooke inequality the specific- 
heat coefficient 7 = c/T appears. For the purpose of dis¬ 
cussing quantum phase transitions it therefore is sensible 
to define critical exponents a and olt that describe the 
critical behavior of the specific-heat coefficient according 
to 7 (r, T = 0) oc |r| _ “, and 7 ( 7 " = 0, T —► 0) <x T~ aT . 
For thermal phase transitions, a obviously coincides with 
the ordinary specific-heat exponent a. For QPTs, we de¬ 
fine a and ar in terms of the susceptibility \r given by 
the second derivative of an appropriate free energy with 
respect to r; the physical meaning of this susceptibility 
depends on the nature of r. The definition of critical 
exponents discussed in this paper is summarized in Ap¬ 
pendix [A] and we will refer to the exponents a, a, /?, 
etc. as the r-exponents, and to ax, or, /?t, etc. as the 
T-exponents. 

The concepts of scaling, and exponent relations, also 
carry over to dynamical critical phenomena!^) A dynam¬ 
ical critical exponent z is defined by how the critical time 



scale T£ scales with the diverging correlation length £: 
77 oc £ z . In the case of thermal phase transitions, the dy¬ 
namics are decoupled from the statics, and the static ex¬ 
ponents can be determined independently of z. At quan¬ 
tum phase transitions this is not true since the statics and 
dynamics are intrinsically coupled. However, as we will 
show, at T = 0 there still are many exponent relations 
that involve the static exponents only. A further compli¬ 
cation is the fact that multiple critical time scales, and 
hence multiple dynamical exponents z, which do occur at 
classical transitional but are not common, are the rule 
rather than the exception at quantum phase transitions. 

Interestingly, the same concepts can also be applied to 
first-order phase transitions. Fisher and Berkci 11 have 
developed a scaling description of classical first-order 
transitions and have shown that they can be understood 
as a limiting case of second-order transitions or critical 
points. 

The purpose of this paper is to thoroughly discuss these 
concepts in the context of quantum phase transitions. 
While scaling concepts have been used since the notion 
of quantum phase transitions was invented, the issue of 
multiple dynamical exponents and its interplay with dan¬ 
gerous irrelevant variables has never been systematically 
discussed. Similarly, the importance of exponent rela¬ 
tions has not been stressed in this context, and even 
the distinction between r-exponents and T-exponents is 
not always clearly made. In the first part of the paper 
we discuss scaling concepts for quantum phase transi¬ 
tions in general, with an emphasis on which aspects can 
be taken over from classical phase-transition theory and 
which ones require modifications or new ideas. In par¬ 
ticular, we show that the specific-heat coefficient shows 
scaling behavior that sets it apart from other observables 
and is not in direct analogy to the classical case, and 
we derive exponent relations that involve the specific- 
heat exponents. Exponent relations have proven very 
useful in the classical context, since they provide strin¬ 
gent checks for whether experiments are actually in the 
asymptotic region. They are expected to be similarly 
useful for QPTs. We also generalized the scaling theory 
for thermal first-order transitions by Fisher and Berkci 11 
to the quantum case. In the second part of the paper we 
apply all of these concepts to the problem of the quan¬ 
tum ferromagnetic transition in metals. This transition 
is known to be generically first order in clean systems, 
and second order in disordered ones, and it has multi¬ 
ple dynamic critical exponents, which allows for an il¬ 
lustration of all of the ideas developed in the first part. 
It also displays many crossover phenomena the under¬ 
standing of which is crucial for the interpretation of ex¬ 
periments. These crossovers can be phrased in terms of 
fixed points, originally discussed by Hertz,® that become 
unstable asymptotically close to the transition. The ex¬ 
perimental relevance of Hertz’s fixed point in disordered 
systems in particular is discussed here for the first time. 

This paper is organized as follows. In Sec. [H] we derive 
exponent relations for the quantum case, both for the 
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usual exponents and for those unique to QPTs, and we 
generalize the notion of scaling at a first-order transition 
to the quantum case. In Sec. |III| we apply our results to 
the problem of the quantum ferromagnetic transition in 
metals, and in Sec. |IV| we conclude with a general discus¬ 
sion. Definitions of common critical exponents are given 
in Appendix [A] and various classical exponent relations 
are recalled in Appendix [B] Appendix [C] summarizes the 
Harris criterion and its generalization for the correlation- 
length exponent in systems with quenched disorder. 


II. SCALING AT QUANTUM PHASE 
TRANSITIONS 

A. General concepts for quantum phase transitions 

In this section we discuss some very general points re¬ 
lated to the role of dangerous irrelevant variables and 
multiple time scales in the context of quantum phase 
transitions. 


1. Dangerous irrelevant variables and dynamical scaling 


The classical exponent relations recalled in Appendix 
[B] fall into two distinct classes: The Widom, Essam- 
Fisher, and Fisher relations hold quite generally, includ¬ 
ing the case of systems above an upper critical dimension¬ 
ality. By contrast, the hyperscaling relation, Eq. (BidI, 


holds only below an upper critical dimension, when no 
dangerous irrelevant variables (DIVs) are present. At 
quantum critical points the role of DIVs is more complex 
than in the classical case. A DIV can affect the temper¬ 
ature scaling of an observable, and it can do so with or 
without also affecting the static scaling. We will refer to 
the assumptions underlying scaling that can or cannot be 
altered by DIVs as “strong” and “weak” scaling assump¬ 
tions, respectively, and to the resulting scaling behavior 
as “weak scaling” and “strong scaling”. We will further 
distinguish between weak and strong static and dynamic 
scaling as appropriate. With this nomenclature, hyper¬ 
scaling in the usual sense is a special case of strong static 
scaling. (Note that weak scaling is more robust than 
strong scaling.) The DIV concept is more important for 
quantum phase transitions than for classical ones, since 
the latter are more likely to be above their upper critical 
dimension and strong scaling is more often violated. 

To illustrate these points, consider an observable O. 
Let [O] be the scale dimension of O as determined by 
power counting within the framework of an appropriate 
field theory, and define a critical exponent oo by the crit¬ 
ical behavior of O as a function of the dimensionless dis¬ 
tance from criticality at T = 0, which we denote by r: 
0(r,T = 0) oc |r|“h If strong static scaling holds, O 
obeys a homogeneity law 12 


0(r,T = 0) = b~{° ] $ 0 (r& 1/l/ ,T = 0) . 


(2.1a) 


Here v is the correlation length exponent, and <f>o is a 
scaling function. Strong scaling thus leads to 

oo = u[0\. (2.1b) 

This does not remain true if there is an irrelevant variable 
that is dangerous with respect to the r-dependence of 
O. This changes the value of oo, but usually it does not 
change the fact that O obeys a homogeneity law. Weak 
static scaling thus still holds in general,^ as expressed 

by 

0(r,T = O) =6- w/, '$ o (r& 1/l ',T = 0) . (2.2) 

The value of oo is now affected by the pertinent DIV and 
can no longer be determined by scaling arguments alone, 
even if [O] and v are known. 

So far we have just stated the usual concept of a 
DIV, EH applied to the static scaling behavior at a quan¬ 
tum phase transition. Now consider the temperature 
dependence of O and define a critical exponent lot by 
0{r = 0, T) oc T UT . Let z be the dynamical exponent 
as obtained from the underlying field theory by power 
counting. (For now we assume there is only one such z, 
we will generalize to the case of multiple dynamical expo¬ 
nents below.) Strong dynamical scaling then is expressed 
by a homogeneity law 

0{r,T) = b~ u / v §o{rb 1 / v ,Tb z ) , (2.3a) 

which yields 

out = oo/vz . (2.3b) 


If an irrelevant variable is dangerous with respect to 
the temperature dependence of O , the relation (2.3b) no 
longer holds. However, as in the static case we still have 
a homogeneity law as long as weak scaling holds. We 
write 


0(r,T) = b~ u / u ^ 0 {rb 1/v ,Tb zo ) , (2.4a) 


where 


zo = u/wt 


(2.4b) 


with oo and out the physical exponents that include the 
effects of the DIV. Equations (2.41 are valid as long as 


weak scaling holds. If in addition strong static and/or 
dynami c scal ing holds one also has the relations (2.1b ) 
and/or (2.3b I between exponents. 

Now consider the free-energy density / as a function 
of r, T, and the held h conjugate to the order parameter. 
Assuming strong scaling, we have 


f{r,h,T) = b~( d+z) <i> f (rb 1/v ,hb [h] ,Tb z ) , (2.5) 


where [h] is the power-counting scale dimension of h. We 
define a “control-parameter susceptibility” Xr as Xr = 
— d 2 f/dr 2 . The physical meaning of depends on the 
nature of the control parameter. For instance, if the QPT 
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is triggered by hydrostatic pressure, yy will be propor¬ 
tional to the compressibility of the system. We further 
define a critical exponent a by y>(r ,T = 0 ) oc |?’| _a . 
Note that the such-defined a at a QPT has nothing to 
do with the specific heat; however, at a thermal transi¬ 
tion our definition makes \ r the specific-heat coefficient, 
and a has its usual meaning. We can now incorporate 
the effect of any DIVs by writing 

f{r,h,T) = b~ i - 2 - a) / 1 ' ^ f {rb l / v ,hb^ /v ,Tb p/vpT ) , 

( 2 . 6 ) 

which is valid as long as weak scaling holds. 


2. Multiple time scales 


So far we have assumed that there is a single underlying 
dynamical exponent z that may or may not be modified 
by a DIV. An additional complication is that often more 
than one dynamical exponent is present. This can hap¬ 
pen at classical phase transitions™ but it is much more 
common at quantum phase transitions, chiefly because 
there are more soft or massless modes at T = 0 than at 
T > 0. For instance, at any quantum phase transition in 
a metallic system the coupling of the order parameter to 
the conduction electrons introduces a second time scale: 
a ballistic one (z = 1 ) in clean systems, or a diffusive one 
(z = 2 ) in disordered onesP^This issue is independent of 
whether or not DIVs invalidate strong scaling, so for sim¬ 
plicity we will assume strong scaling in this subsection. 

For definiteness, let as assume that there are two criti¬ 
cal time scales with dynamical exponents z\ and Zi < Z\. 
The generalization of the homogeneity law (2.1a) for O 
then reads 

0(r,T) = b~ [ ° ] § 0 {rb d/v ,Tb zl ,Tb Z2 ) . (2.7) 

Let us assume [O] > 0 (for [0\ < 0 obvious modifica¬ 
tions of the following statements hold). Then the lead¬ 
ing temperature dependence of O at r = 0 will be given 
by z±, i.e., u>t = \0\/zi. However, there is no guaran¬ 
tee that <I>(0 ,:e,?/) depends on x. If it does not, then 
wt = \ 0 ]/z 2 , and the singular T-dependence of O at 
criticality is weaker than one would naively expect. We 
will see explicit examples of this in Sec. |III| 

For the free energy, multiple dynamical exponents re¬ 
sult in multiple scaling parts. Let us again consider the 
case of two critical time scales with dynamical exponents 
z\ > Z 2 - All of the following considerations can be triv¬ 
ially simplified to the case z\ = z%, and they are easily 
generalized to the case of more than two dynamical ex¬ 
ponents. We generalize Eq. (2.51 to 


f(r,h,T) = b-^ 1 ' 1 ^f\rb 1/,y ,hb [h] ,Tb Zl ,Tb Z2 ) 

+ b~(d+Z2) $(2) ( r b l/^ h b [h ] ) T b Z2 ) _ 

( 2 . 8 ) 

( 2 ) 

Note that the scaling function ’ does not depend on 
the argument T b Zl . This is because the entropy den¬ 


sity, s = df/dT , must vanish at least as fast as l/£ d 
for t; —> ocP-2 For / and its derivatives with respect to r 
and h the second term obviously yields the most singular 
contribution. We thus can obtain homogeneity laws for 
the order parameter m = df/dh, the order-parameter 
susceptibility Xm = d 2 f/dh 2 , and the control-parameter 
susceptibility \r = d 2 f /dr 2 from a scaling part of the 
free-energy density 

f 2 (r,h,T) =b- {d+Z2) $ { f\rb 1/v ,hb [h \Tb Z2 ) . (2.9) 

For the specific-heat coefficient 7 = c/T = —d 2 f/dT 2 , 
on the other hand, the largest z yields the most singular 
contribution and we have 


7 (r, T) =b Zl ~ d <f> 7 (r b 1/v , T b Zl ,Tb Z2 ) . 


( 2 . 10 ) 


The leading temperature dependence at r = 0 is in gen¬ 
eral given by the first argument of the scaling function, 
i.e., y(r = 0,T) oc T d ^ Zl ~ 1 . 

We emphasize again that all of the above relations as¬ 
sume strong scaling. See below for a discussion of the 
weaker statements that remain valid if strong scaling is 
violated. We also stress that Eq. (2.8), with its two addi¬ 


tive scaling parts, is an ansatz. In the case of the quan¬ 
tum ferromagnetic transition discussed as our prime ex¬ 
ample in Sec. |III| it is known explicitly that the scaling 
part of the free energy has this form. 


B. Exponent relations at quantum critical points 

We now use the concepts laid out above to derive var¬ 
ious relations between critical exponents. 


1. Exponent relations that rely on weak scaling 


Let the observable in question be the order parameter 
to, let h be the field conjugate to the order parameter, 
and let x = dm/dh be the order-parameter susceptibility. 
As long as weak scaling holds, Eqs. (2.41, generalized to 
allow for the dependence of to on h, will be valid with 
u> = /? and lot = Pt- 

m{r,h,T) = b~ ll> / lJ ^ m {rb llv ,hb p5 / v ,Tb Zm ) , (2.11) 


where z m = P/v^t is the dynamical exponent relevant 
for the order parameter. In the context of Eq. (2.8 1 , z m 
is Z 2 , possibly modified by DIVs. Differentiating with 
respect to h we find a corresponding homogeneity law 
for the order-parameter susceptibility: 

Xm (r,T) = ( r 6 1 /*',r 6 *“) . (2.12) 


With the definitions of the exponents 7 and 7 ^, Appendix 
[AJ this implies the 

Widom equality: 

7 = p(8-l) , 

7 t = Pt{5 — 1 ) , 


(2.13a) 

(2.13b) 
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which holds at a QPT for both the r-exponents and the 
T-exponents as a consequence of weak scaling only. It is 
not affected by the presence of multiple time scales. 

Next we consider, in addition to to, the order- 
parameter susceptibility Xm and the control-parameter 
susceptibility y r , which all are obtained from the scaling 
part f -2 of the free-energy density, Eq. Incorporat¬ 

ing the effects of DIVs, if any, as in Sec. IIA 1[ we can 
write 


m(r , h, T) 


2)1 v <& m {rb 1/v ,hbP slv ,T 

b Zm ) , 




(2.14a) 

Xm(r, T) 

_ ^(2/35 +a 

- 2)/!/ $ x .(t 'b 1 / v ,Tb z ™) , 

(2.14b) 

Xr(r, T) 

= b a/u A> r 

(rb 1/v ,Tb z ™) . 

(2.14c) 

From the definitions of 

the exponents (3, 7, (3t, 

7 t, and 


or, and using Eq. (2.13a I we then obtain the 


Essam-Fisher equality: 

a + 2(3 + 7 = 2 , 

for the r-exponents, and its analog 

Ut + 2 Pt + 7 t = 2 /vzr, 


(2.15a) 


(2.15b) 


for the T-exponents. Equation ( |2. 15a ) depends on weak 
scaling only. Note that a is not the specific-heat ex¬ 
ponent, but rather the control-parameter-susceptibility 
exponent defined in Sec. |II A 1| and Appendix [A] Equa¬ 
tion (2.15a) is the natural extension of the classical 
Essam-Fisher equality to quantum phase transitions, 
even though it does not involve a specific-heat exponent. 
The simple relation between Eq. (2.15aI and (2.15bI, 


which is consistent with Eq. (2.4b I with zq = z m , is due 


to the fact that the dominant dynamical exponent is the 
same for all three observables to, Xmi an d Xr- The latter 
can be guaranteed only in the presence of strong scaling, 
see the remarks after Eq. ( |2.8|). However, while strong 
scaling is su fficient for Eq. ( |2.15b I, it is not necessary. 

Equation (2.14a) also implies — (3 = (35 + a — 2. To¬ 


gether with the Essam-Fisher equality this allows us to 
express 5 in terms of a and 7: 


<5 = 


2 — a 


.7 


2 — a — 7 


(2.16) 


This extends to QPTs another relation that is well known 
for thermal phase transitions!- 1 Note that Eq. (2.161 is not 


independent; it follows from a combination of the Widom 
and Essam-Fisher equalities. We also stress again that in 
the current context a is not the specific-heat coefficient. 

Now consider the order-parameter two-point correla¬ 
tion function G(x) = dr (m(x,r) m( 0,0)), with r 
the imaginary-time variable. The spatial integral of 
G{x) yields the order-parameter susceptibility, Xm = 
f dx G(x). It has the form 


G(x) = 


o- M/5 

r irf—2+77 1 


(2.17) 


which defines both the correlation length £ and the crit¬ 
ical exponent 77. One thus has 


Xm °c £ 


2 —v 


(2.18) 


An equivalent argument is to generalize the homogeneity 
equation for the order-parameter susceptibility to include 
the wave-number dependence: 

Xm (r,T-k) =b^^ x (rb 1 /’',Tb z -,kb) . (2.19) 

With the definitions of the exponents v and vt, Ap¬ 
pendix [A] we obtain from either Eq. (2.18) or (2.191 the 


Fisher equality: 


7 = (2 — r))v , 
7t = (2 — 77) 


(2.20a) 

(2.20b) 


for both the r-exponents and the T-exponents. It de¬ 
pends on weak scaling only. 

In summary, we have found three independent weak- 
scaling exponent relations for the r-exponents, viz. the 
Widom, Fisher, and Essam-Fisher equalities (2.13aI, 


(2.15aI, and (2.20aI, and three corresponding relations 


for the T-exponents, Eqs. (2.13b), (2.15bI, and (2.20bI. 


2. A rigorous inequality 

Exponent relations that involve the specific-heat ex¬ 
ponent a (see the Introduction and Appendix [A] for a 
definition and discussion of a) fall into a different class, 
since the critical behavior of the specific-heat coefficient 
is in general governed by a scaling part of the free en¬ 
ergy that is different from the one that determines the 
exponents discussed so far, see Eq. ( 2. 10| ) and the related 
discussion. We start with the thermodynamic identity 
that underlies the classical Rushbrooke inequalitjffG^ 


1 - Im/lh = [{dm/dT)l] /xTlh ■ 


( 2 . 21 ) 


Here and denote the specific-heat coefficient at 
fixed order parameter and fixed conjugate field, respec¬ 
tively, and xt is the isothermal OP susceptibility. Note 
that, in order to apply this relation to QPTs, it is crucial 
to formulate it in terms of the specific-heat coefficients; 
the usual formulation in terms of the specific heats leads 
to a factor of T on the right-hand side that makes the 
T — > 0 limit ill defined. 1 '^ With the T-exponents as de¬ 
fined in Appendix [A] this yields an exponent inequality 
at a QPT in the form 


OiT + 2/3t ~t~ 7 t (3 2 . 


(2.22a) 


This is rigorous, since it depends only on thermodynamic 
stability arguments. If in addition we use weak scaling of 
the order parameter, Eq. (2.111, we find a corresponding 


inequality for the 7"-exponents, 

a + 2/3 + 7 > 2(3/[3t ■ 


(2.22b) 
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This follows since weak scaling of the OP implies 
( dm/dT)T-o,h~o oc (—r)^ 1-1 /^. Note that this is 
different from the classical Rushbrooke inequality, Eq. 
( |B5| ). As for the latter, Eqs. (2.221 hold as equalities 
if lm/"1h —> 1 for T —> 0 at r = 0, and for r —► 0 at 
T = 0, respectively, and as inequalities otherwise. In 
the classical limit, both of the Eqs. (2.22) turn into the 
the classical Rushbrooke inequality: The T-exponents co¬ 
incide with the r-exponents, and a coincides with the 
classical specific-heat exponent a, see the discussion in 
Appendix [A] 


3. Hyperscaling relations 


The exponent relations we derived and discussed so far 
depended at most on weak scaling (with the exception 
of Eq. (2.15b), which relies to some extent on a strong¬ 


scaling assumption). If strong scaling is valid we can 
derive additional constraints on the exponents. For this 


purpose, we return to Eq. (2.81. As we discussed in this 


context, Z 2 = z m is the dynamical exponent for the or¬ 
der parameter, the order-parameter susceptibility, and 
the control-parameter susceptibility, while z\ = z c > z m 
is the dynamical exponent for the specific-heat coeffi¬ 
cient. In general z c could be z\ modified by DIVs, but 
the strong-scaling assumption means z c = Z\. Strong 
scaling thus implies the quantum versions of the usual 
hyperscaling relation (BidI, 

(2.23a) 
(2.23b) 


a = 2 - v{d + Zm) , 
a T = 2,/vz m - d/z m - 1 . 


The latter in agreement with Eq. (2.3bI, and the former 


(for the case where there is only one dynamical criti¬ 
cal exponent) with Ref. (T5] We will refer to these as 
the a-hyperscaling relations. Analogously, we find from 
the first term in Eq. (2.81 hyperscaling relations for the 
specific-heat exponents 


(2.24a) 
(2.24b) 

which we will refer to as the d-hyperscaling relations. 
Equations (2.231 and (2.241 are two independent 


a = is(z c - d) , 
olt = 1 - d/z c , 


strong-scaling relations. In conjunction with the weak- 
scaling relations from Sec. |IIB l| we can use them to de¬ 
rive additional relations that are not independent. For 


instance, Eqs. (2.231 and (2.241 imply a relation between 
the control-parameter exponent a and the specific-heat 
exponent d, 

a = a — 2 + i/(z m + z c ) , (2.25a) 

ap = olt - 2/vZm + 2 + d(l/z m - l/z c ) . (2.25b) 


Combining these with the Essam-Fisher relations (2.15) 
we find 

a + 2/3 + 7 = v(z m + z c ) , (2.26a) 

olt + 2 / 3 t + 7 t = 2 + d(l/z m - l/z c ) , 

(2.26b) 


Both of these equalities are consistent with Eqs. (2.221, 
since z c > z m always. Note that they are physically 
different from the Essam-Fisher equalities (2.151, which 
depend on weak scaling only. 

We can use Eq. (2.25a I together with Eq. (2.16) to 
express 5 in terms of d, 7, and the dynamical exponents, 


v(z m + z c )-a + 'y 
v(z m + z c ) - a - 7 ' 


(2.27) 


We note that all of the above relations rely on strong 
scaling, even if they do not explicitly involve the dimen¬ 
sionality. 

Finally, the electrical conductivity is dimensionally an 
inverse length to the power {d — 2). A strong-scaling 
hypothesis thus implies that the scaling part Act of the 
conductivity obeys a homogeneity law 

Acr(r, T) = &-(<*- 2 ) F a (rb 1/v ,Tb z c) . (2.28) 


For the exponents s and St defined in Eq. (A6) this im¬ 
plies the 


Wegner equality: 


s = v(d — 2) , (2.29a) 

s T = (d-2)/z c . (2.29b) 

These relations play an important role in the theory of 
electron localization) 19 20 but are much more generally 
applicable. They do, however, depend on strong scaling 
and are not valid if DIVs affect the conductivity. Note 
that in general, with DIVs taken into account, s and 
St can be either positive or negative in any dimension, 
with negative values applying to either the clean limit or 
regimes where the residual resistivity is small compared 
to the temperature-dependent part. 21 


C. Scaling at quantum first-order transitions 

Fisher and Berker have shown how a classical first- 
order transition can be understood within a standard 
scaling and RG framework. 1 ! In a RG context, at a first- 
order transition all gradient-free operators in a Landau- 
Ginzburg-Wilson (LGW) theory are relevant and have 
the largest scale dimension that is thermodynamically 
allowed, viz., d. This just follows from the fact that the 
OP itself is dimensionless at a first-order transition, so 
the operator dimensions must make up for the scale di¬ 
mension of the spatial integration measure i 22 -! In partic¬ 
ular, v = 1/d. This is reconciled with the usual notion 
of a finite correlation length at a first-order transition by 
means of finite-size scaling considerations. A dimension¬ 
less order parameter requires [h] = d, which is equivalent 
to /3 = 0, and the remaining classical exponent values 
follow readily: 77 = 2 — d, 7=1, S = 00 , and a = 1. All 
scaling relations, Appendix [5] are fulfilled, including the 
hyperscaling relation. (Note that DIVs are not an issue 
in this case, since the suspect operators are relevant.) 
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The reasonin g of Fisher and Berker can readily be ex¬ 
tended to QPTs.GEO For the case of making the logic clear, 
we first do so for a QPT with only one dynamical expo¬ 
nent, and then generalize to the case of two dynamical 
exponents. Consider again the homogeneity law for the 
free energy, Eq. (2.5). In order for the order parameter 
to = df/dh to be dimensionless we must have 


[h]=d + z. (2.30) 

A dimensionless OP in turn implies 

\/v = d + z , (2.31a) 

which generalizes the classical 1/zz = It further im¬ 
plies 


77 = 2 — d — z , S = 00 , /? = 0 , Pt = 0 . 

(2.31b) 

Considering the OP susceptibility Xm = d 2 f/dh 2 , and 
using Eqs. (2.301 and (2.31aI, we have 


7 = 1 


7 t = {d + z)/z . 


(2.31c) 


zero-temperature limit, where the size in the imaginary- 
time direction becomes infinitely large. In this case argu¬ 
ments identical to those used in classical finite-size scaling 
theorjES lead to a time scale that scales as 

t ~ e aLd , (2.32) 

with a a positive constant. In the classical case the analo¬ 
gous result has led to a number of remarkable conclusions 
that have been confirmed experimentally.^ The physical 
meaning of r in the quantum case is the time it takes 
for a droplet of volume L d containing a certain state to 
transform into a different state via a tunneling process. 

Now we turn to the case of two different dynamical 
exponents for the OP and the specific heat, which we 
again denote by z m and z c , respectively. The imaginary- 
time integral in the mass term and the Zeeman term in 
a LGW functional will then have a scale dimension of 
— z m (taking the scale dimension of a length to be — 1 ). 
Making the OP dimensionless thus requires 

[h] = d + z m , (2.33a) 


Similarly, we find from the homogeneity law for the 
control-parameter susceptibility Xr = d 2 f/dr 2 the ex¬ 
ponents 


and 

1 /v = d + z m . (2.33b) 


a = l , otT = (d + z)/z. (2.31d) 

The value a = 1 implies that the first derivative of the 
free energy with respect to r is discontinuous, and the 
second derivative contains a delta-function. This is anal¬ 
ogous to the entropy being discontinuous, and the exis¬ 
tence of a latent heat, at a thermal phase transition. In 
the quantum case the physical interpretation of a = 1 
depends on the nature of the control parameter; see the 
discussion in Sec. |IV B 1| From the homogeneity law for 
the specific-heat coefficient 7 = d 2 f /dT 2 we find 


The OP is now dimensionless by construction, which im¬ 
plies 

rj = 2 — d — z m , 5 = 00 , /3 = 0 , Pt = 0 . 

(2.33c) 

From the OP susceptibility, Xm = dm/dh , we obtain 

7 = 1 , 7 T = (d +z m )/z m , (2.33d) 

and from the control-parameter susceptibility Xr = 
d 2 f/dr 2 


z — d 

z + d 


z — d 

olt = - 

z 


(2.31e) 


Finally, Eq. (2.2 1 holds in particular for the correlation- 
length exponent, which yields 


v T = l/z . (2.31f) 

We see that all scaling relations discussed in Sec. m 
hold, in analogy to the classical case. 

In the above discussion we have assumed what is called 
a block geometry in the finite-size scaling theory of classi¬ 
cal first order phase transitions. Physically it is realized 
in a finite system whose linear sizes in all dimensions, in¬ 
cluding the imaginary-time or inverse temperature one, 
are comparable. Another choice is a cylinder geometry 
in which the linear size in one dimension is large com¬ 
pared to the others. In the quantum case this is of spe¬ 
cial interest for a system with fixed finite volume in the 


a = 1 


a T = (d + z m )/z n 


(2.33e) 


Since the correlation length is entirely a property of the 
OP-OP correlation function, its temperature dependence 
is also governed by z m and Eq. (2.4b I holds for ui = v 
with z = z m : 


v T = 1 /z m ■ (2.33f) 

The specific-heat coefficient is governed by the dynamical 
exponent z Cl and we thus have 

a = is(z c - d) = ~ c ^ , d r = 1 -d/z c . (2.33g) 

All of the scaling relations discussed in Sec. |II A| still hold. 

We will see explicit examples for the general relations 
developed here in Sec. |III| and we will further discuss the 
underlying concepts and their consequences in Sec. m 
















III. APPLICATION TO QUANTUM 
FERROMAGNETS 

The scaling theory developed in Sec. |TT] is in general 
quite complex because of the occurrence of multiple time 
scales at many quantum phase transitions. As an explicit 
illustration we now apply the theory to the quantum fer¬ 
romagnetic phase transitions in metallic systems. For 
a review of the ferromagnetic quantum phase transition 
problem, see Ref. f271 


A. Clean systems 

The ferromagnetic transition in metals was one of the 
earliest quantum phase transitions considered, and it 
was used as an example by Hertz in his seminal paper 
on a renormalization-group approach to quantum phase 
transitions P Hertz considered a ferromagnetic order pa¬ 
rameter m coupled to conduction electrons. The latter 
lead to Landau damping of the magnetic fluctuations, 
and on general grounds, using symmetry arguments as 
well as the known structure of the Fermi-liquid dynam¬ 
ics, one can write down an action 

A H ertz = - y~] [r + ak 2 + c|Q n |/|fc|] 
k,n 

xm(k , Q n ) ■ m(k , —f l n ) + 0(m 4 ) , (3.1) 

where we have written only the term quadratic in the or¬ 
der parameter explicitly. Here fl n is a bosonic Matsubara 
frequency, and r, a, and c are parameters of the LGW 
functional. The term |Q„|/|fe| reflects the dynamics of 
the conduction electrons. This is the action considered 
by Hertz, who explicitly derived it from a specific model. 
From this he concluded that the dynamical critical ex¬ 
ponent is z = 3 (since H ~ |fc| 3 for r = 0), and that the 
upper critical dimension is =4 — 2=1. This in turn 
led to the conclusion that the quantum phase transition 
is second order and described by a simple Gaussian fixed 
point, and that the static critical behavior is mean-field- 
like for all d > 1 with logarithmic corrections to scaling 
for d = 1.® The behavior at finite temperature was con¬ 
sidered in detail by Millis,® and the results of the RG 
treatment for d = 3 confirmed results obtained earlier by 
Moriya and co-workers by means of what is often called 
self-consistent spin-fluctuation theory.® This combined 
body of work is often referred to as Hertz-Millis-Moriya 
theory. 

It was later shown that the above conclusions do not 
hold due to properties of the conduction electrons that 
are not captured in Hertz’s action. A careful analysis of 
the soft fermionic modes coupling to the magnetization 
shows that, in three-dimensional systems, the leading 
wave-vector dependence in the Gaussian action is a term 
proportional to fc 2 ln(l/|fc|) with a negative prefactor® 
and such a term is indeed generated under renormal¬ 
ization if one starts with Hertz’s action. An equivalent 


statement is that there is a to 4 ln(l/m) term in a general¬ 
ized Landau theory for the ferromagnetic quantum phase 
transition. 2 ' As a result, Hertz’s fixed point is not stable, 
and the ferromagnetic quantum pha se transition in clean 
metals is generically first order! 32 33 However, depending 
on quantitative details involving the strength of the elec¬ 
tron correlations and the inevitable weak disorder, there 
can be a sizable regime where the critical behavior asso¬ 
ciated with Hertz’s fixed point is observable, even though 
asymptotically close to the transition it crosses over to 
a first-order transition.® In this subsection we therefore 
discuss, and include in Table |T] the critical exponents 
at Hertz’s fixed point, together with the exponents as¬ 
sociated with the ultimate first-order transition, and the 
critical behavior at the termination point of the tricritical 
wings that result from the first-order transition.® 


1. Hertz’s fixed point 


Table |T] lists the critical exponents associated with 
Hertz’s fixed point. The values for the static expo¬ 
nents are straightforward generalizations of the results 
obtained in Refs. H2S The value for the resistivity ex¬ 
ponent st is Mathon’s result® generalized to d dimen¬ 
sions, and z c = 3 is Hertz’s value for the dynamical 
critical exponent P There was no notion of a separate 
dynamical exponent z m in the original work. However, 
this concept naturally arises if we take the point of view 
expressed in Eqs. (2.41 and incorporate the effect s of 
the DIVs in the homogeneity laws. Equation ( |2.4b I ap¬ 
plied to the order parameter then implies z m = /3/iz/St, 
which yields z m = 1//3t as listed in Table [I] Alterna¬ 
tively, one can work with only one dynamical exponent, 
2 = 3, and consider the DIV explicitly, as was done by 
Millis,®see also Ref. 1271 This reasoning effectively leads 
to z m = 2/(1 + v{d — 1)), which is the same result as 
above. At the upper critical dimension d = 1 the effects 
of the DIVs disappear, and there is only one dynami¬ 
cal exponent z m = z c = 3. Note that there is only one 
fundamental time scale in Hertz’s theory, and the ap¬ 
pearance of a second dynamical exponent is solely due 
to the fact that the fundamental 2 is modified by a DIV 
in some contexts, but not in others. This is very differ¬ 
ent from the theories discussed in Secs. |III A2| |IIIA3| 
and |IIIB 2[ which intrinsically contain two time scales 
of different physical origin. This is an indication that 
important physics related to the conduction electrons is 
missing in Hertz’s theory, see Sec. |IV B 2| for additional 
comments on this issue. 

Since for d > 1 the system is above its upper critical di¬ 
mension, DIVs in general invalidate any scaling relations 
that rely on strong scaling as defined in Sec. |II A| while 
those that rely on weak scaling only will remain valid. 
Indeed, considering the values in the table, we see that 
the Widom equalities (2.131, the Essam-Fisher equalities 


(2.151, the Fisher equalities (2.201, and the Rushbrooke 


inequalities (2.22) are satisfied for all d > 1. The hyper- 





























9 





Figure 1: Hertz-type scaling behavior as observed in Ni3Ali_ :E Ga2. From left to right: T c vs x phase diagram, inverse magnetic 
susceptibility vs. T 4 ^ 3 for the critical sample, magnetization squared vs. T 4 ^ 3 for various concentrations. Figure adapted from 
Ref. En¬ 


sealing relations (2.23), on the other hand, do not hold. 


The hyperscaling relations (2.241 are a more complicated 
case. The values of a and cxt listed in Table |T] reflect the 
leading fluctuation contribution to the specific-heat coef¬ 
ficient, which do obey strong scaling. For 1 < d < 3 they 
dominate the constant mean-field contribution which vi¬ 


olates strong scaling, and as a result Eqs. (2.24) hold even 


though the system is above its upper critical dimension. 
For d > 3 they are subleading compared to the mean-field 
contribution, the specific-heat exponents lock into their 
mean-field values a = olt = 0, and hyperscaling breaks 
down. 

The contribution to the resistivity p due to scatter¬ 
ing of electrons by critical fluctuations was calculated by 
MathonP who found A p(r = 0,T) oc T 5 / 3 in d = 3. A 
simple generalization to general d yields A p(r = 0, T) oc 
rp(d+ 2)/3, or St = — (^ _|_ 2)/3 for the conductivity expo¬ 
nent St defined in Appendix [A] From a scaling point 
of view, this result can be made plausible as follows. 


Consider the strong-scaling homogeneity law (2.28) for 


the conductivity contribution Act. In a clean system, 
the backscattering factor in the Boltzmann equation pro¬ 
vides an additional factor of the hydrodynamic momen¬ 
tum squared, which is not captured by power counting. 
This leads to a scale dimension of Act that is effectively 
equal to d — 4 rather than d — 2. In addition, the DIV 
u with scale dimension [u] = — (d — 1) affects Act, and 
from Fermi’s golden rule it is plausible that Act oc 1/it 2 . 
Combining these arguments, we have 

Aa{r,T) = &- (d - 4) 6 2(d " 1} F CT (r 6 1/!/ , T & Zc ) 

= & (d+2) F <J (rb 1/v ,Tb z °) , (3.2) 

which yields st = ~{d + 2)/z c as quoted above and listed 
in Table [I] Note that the Wegner equality (2.29b) is vio¬ 
lated for two reasons: First, the conductivity, or the un¬ 
derlying relaxation rate, does not have its power-counting 
scale dimension due to the backscattering factor in the 
Boltzmann equation. Second, the DIV further modifies 
the scale dimension and actually changes its sign for all 
dimensions. Finally, note that the exponent s as defined 


in Eq. (A6l does not exist since a(T = 0) = oo. Rather, 


for T —> 0 at fixed r yf 0 the conductivity will cross over 
to the Fermi-liquid result, Act(T) oc T~ 2 , and Eq. (3.2 1 
yields the r-dependence of the prefactor: 

A<r(r yf 0, T —>■ 0) oc |r|( 4-d)/2 T~ 2 . (3.3) 


As an illustration, we discuss an experiment on the 
ferromagnet NRAR^Ga,,;, which displays a QPT at 
x = x c « 0.34p3 This system is only moderately dis¬ 
ordered, and critical behavior associated with the clean 
Hertz fixed point is expected to be observable in a siz¬ 
able transient regime. 27 There are three experimental ob¬ 
servations, from which four exponents can be deduced: 

(1) The critical temperature scales as T c ~ r 3 / 4 , see 
the first panel in Fig. [l] This implies vz m = v z/{ 1 + 
2v) = 3/4, as can be seen, for instance, from Eq. (2.121: 


The critical temperature is determined by Xm diverging, 
which must happen for a particular value of the argument 
x in $ x (l,cc). This in turn implies T c /r UZm = const. 

(2) At the critical concentration, the magnetic suscep¬ 
tibility scales as Xm ~ T -4 / 3 , see seconc j panel in 
Fig. 0 This implies 7 t = 4/3. Together with the first 
observation, it also implies 7=1, since r ~ T 4//3 ~ 1 /\. 

(3) The magnetization vanishes as m 2 oc T 4 ^ 3 — T 4 / 3 , see 
the third panel in Fig. [T] In addition to confirming the 
product vz m = 3/4 this yields /3 = 1/2, see the discussion 
after Eq. (4.21. All of these results are in agreement with 
the theoretical results summarized in Table |T] specialized 
to d = 3. The conductivity exponent st = —5/3 has also 
been observed in various materials, see, e.g., Ref. 00] 


2. The first-order transition 

There are strong theoretical arguments for the quan¬ 
tum phase transition from a param agnet to a homoge¬ 
neous ferromagnet to be first order OSES] anc j jqqg j s j n _ 
deed the prevalent experimental observation. 2 ' Here we 
discuss the exponent values, and the scaling relations, at 
this particular first-order quantum phase transition as an 
example of the general scaling theory in Sec. |II C] 
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Table I: Static, dynamic, and transport critical exponents as defined in the text, and validity of exponent relations, at various 
fixed points for the quantum phase transition in metallic ferromagnets. The first two columns of values represent true asymptotic 
exponents (QWCP = quantum wing-critical point). The third and fourth columns of values refer to critical behavior that is 
asymptotically not represented by pure power laws, and the fifth and sixth columns of values represent power-law behavior 
that is observable only in a transient non-asymptotic regime. Strong-scaling relations for the T-exponents are not shown, they 
hold if an only if the corresponding relations for the r-exponents hold. See the text for additional explanation. N/A = not 
applicable. 
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0 
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O 
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a 

1 

0 
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0 

0 

x 
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“For d > 1. 

b This physical fixed point maps onto the unphysical (describing pre-asymptotic behavior only) clean Hertz fixed point. See 
the text for a discussion of what is observed if the critical point is approached along generic paths in the phase diagram. 
“For d > 1. At the upper critical dimension d = 1 there are logarithmic corrections to scaling. 

d For 2 < d < 4. Values in this column equal values in the next column for A = 0. The critical behavior consists of log¬ 
normal terms multiplying power laws with the exponents shown. 

“ For 2 < d < 4. A depends on the distance from criticality. In d = 3, A ~ 2/3 in a large region. 

1 For d > 0. At the upper critical dimension d = 0 there are logarithmic corrections to scaling. 


9 For 1 < d < 3 a nd d > 3, respectively. For 0 < d < 4 and d > 4, respectively. 

for the sense in which this result is valid. 9 Except for d = 1. 


i 


IIIB1 


See Sec. 

For 1 < d < 3 only. m For 0 < d < 4 only. n Not an independent relation. 


Except for d = 0. 

See the text for a discussion. 
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Let us first discuss the dynamical exponents. The soft 
fermionic fluctuations that drive the transition first order 
are of a ballistic nature with z = z% = 1. Their coupling 
to the order-parameter fluctuations lead to z m = 1. The 
static fermionic susceptibility in clean metals at T = 0 
has a wave-number dependence y(fc —► 0) oc const. + 
|fc| rf_ 1 EH This plays against the |fi|/|fc| Landau-damping 
term (see Eq. (3.11 with ak 2 replaced by |fc| d_1 ), which 
produces another dynamical exponent Zi = d. For d > 1 
we thus have z c = Z\ = d. 0 = /3t = 0 and S = 00 
by the discontinuous nature of the transition. Also, the 
fact that the order parameter is dimensionless enforces 
v = 1 /(d + Zm) = l/( d+ 1 ) and r) = 2 — d — z m = 1 — d, 
as explained in Sec. II C[ while ut = vfvzm = 1. It 
further implies that the free-energy density scales linearly 
with the control parameter, which implies a = 1 and 
a.T = a/vz m = d + 1. Finally, in a fermionic system 
the specific-heat coefficient must display a discontinuity 
at a discontinuous phase transition, which implies a = 
ut = 0. Note that, combined with the scaling result 
(2.33g), this implies that the dynamical exponent z c must 
be equal to d. 

All of these exponent values are displayed in Table [T] 
We see that all exponent relations, including those that 
rely on strong scaling, are valid as expected. The Widom 
equation is obviously fulfilled only in the sense that 0(6 — 
1) i 1 as 0 —^ 0, S —y 00 . 


3. The quantum wing-critical point 

Since the quantum ferromagnetic transition is first or¬ 
der, while the corresponding thermal transition is gener- 
ically second order, there necessarily is a tricritical point 
in the phase diagram if the Curie temperature is contin¬ 
uously suppressed by means of some control parameter. 
This in turn leads to the existence of tricritical wings, i.e., 
a pair of surfaces of first-order transitions that emanate 
from the coexistence curve in the h = 0 plane PH These 
wings are indeed commonly observed, see Fig. [2] for an 
example. They end in a pair of quantum critical points 
in the T = 0 plane. These quantum wing-critical points 
(QWCPs; in the literature they are often erroneously re¬ 
ferred to as quantum critical end points) correspond to a 
fixed point that maps onto Hertz’s fixed point; they thus 
are an example of a true quantum critical point that is 
correctly described by Hertz theory, and the critical be¬ 
havior is known exactly for all d > 1 GUI The critical expo¬ 
nents are thus the same as those discussed in Sec. |III A~T| 
and the entries in Table [I] reflect this. However, the appli¬ 
cation of these results for the prediction of experimental 
observations must be handled with care, as we will now 
discuss. 

A crucial point to remember is that at the QWCP, in 
contrast to Hertz’s fixed point, the field h conjugate to 
the order parameter is not the physical magnetic field, 
and h = 0 only on special paths in the phase diagram 
that are not natural paths to choose for an experiment. 


Let H be the physical magnetic field, and for definiteness 
let us assume that the control parameter is hydrostatic 
pressure p. In the three-dimensional T-p-H phase dia¬ 
gram, let the QWCP be located at (0,H c ,p c ), and let 
the magnetization at this point have the value m c . Then 
the conjugate field is given by h = 2 m c 6p — SH , where 
6p = p — p c , and SH = H — jjf c J 27 l 34 l I n order to observe 
the exponent 0, for instance, one therefore needs to ap¬ 
proach the QWCP on a curve given asymptotically by 
SH = 2m c Sp, see path (1) in Fig. 2j This is in exact anal¬ 
ogy to the case of a classical liquid-gas critical point in the 
p-T plane, where an observation of 0 requires that the 
critical point be approached on the critical isochore 
Probing the QWCP along a generic path in the T = 0 
plane measures the exponent <5 instead, e.g., 

m(p = p c ,H,T = 0) oc \H-H c \ l,s = \H-H C \ 1/3 . (3.4) 

This is the behavior predicted for the magnetization 
along path (2) in Fig. [2] A related complication oc¬ 
curs for the temperature dependence of observables at 
the QWCP. Consider again the order parameter. The 
exponent 0t is defined via the T-dependence of m at crit¬ 
icality in zero conjugate field, see Appendix [Aj Since p c is 
temperature dependent, \p c (T) — p c (T = 0)| a T ( d + 1 )/ 3 , 
see Sec. |III A 1| this means that 0t can be observed only 
on a particular surface in the three-dimensional parame¬ 
ter space spanned by T, p, and H. Along a generic path 
through the QWCP the temperature dependence of m is 
given by the exponent combination 2 u T /S. In particu¬ 
lar, just raising the temperature from zero at the critical 
point (path (3) in Fig. [2j yields 

m(p = p c , H = H c , T) oc \h\ 1/s = -T 2lAr / 5 = _T (d+1)/9 . 

(3.5a) 

This is the result that was derived in Ref. [3H see also 
Ref. 23 An observable that is easier to measure is the 
magnetic susceptibility \m , which along the same path 
behaves as 

Xm(p =Pc,H = H c , T) cx I h\-^P s = 7 1 “ 2 ( d + 1 )/ 9 . 

(3.5b) 


B. Disordered systems 

Quenched disorder in metallic ferromagnets introduces 
many different effects. Some, such as the change of the 
conduction-electron dynamics from ballistic to diffusive, 
directly affect the behavior at the quantum phase transi¬ 
tion. Others, such as rare-region effects that can lead to 
the appearance of a q uantum Griffiths region in the para¬ 
magnetic phasep2l32 are superimposed on critical singu¬ 
larities and may easily be confused with the latter. Dis¬ 
entangling these various effects is challenging from both a 
theoretical and an experimental point of view, and many 
open questions remain. Here we ignore rare-region effects 
and discuss the effects of quenched disorder on the phase 
transition itself. 
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Figure 2: Tricritical wings observed in the temperature- 
pressure-magnetic field phase diagram of UGe 2 by Kotegawa 
et al., Ref. El Data points are represented by squares, blue 
planes are planes of first-order transitions, red solid lines are 
lines of second-order transitions, red points are the quantum 
wing-critical points (QWCPs). Also shown is the tricritical 
point (TCP). Dotted green arrows labeled (1), (2), (3) are 
three different paths through a QWCP; the critical behav¬ 
ior predicted for each path is discussed in the text. Figure 
adapted from Ref. El 


This problem was also considered by Hertz, El who ar¬ 
gued that the only salient change compared to the clean 
case is in the Landau-damping term in Eq. (3.1), which 
now is \n n \/k 2 due to the diffusive dynamics of the con¬ 
duction electrons. This obviously leads to a dynamical 
critical exponent z = 4, and to an upper critical dimen¬ 
sion d[ j~ = 0. The finite-temperature behavior can be 
discussed in exact analogy to Millis’s treatment of the 
clean case in Ref. [53 


Hertz’s fixed point in the disordered case is again un¬ 
stable. The physics behind this instability is the same 
as in the clean case, viz. the coupling of the magne¬ 
tization to fermionic soft modes, which now are diffu¬ 
sive in nature at asymptotically small wave numbers. 
At larger wave numbers they cross over to the clean 
soft modes. For weak disorder, this crossover occurs 
at a very small wave number, and one expects the ef¬ 
fects to be small. This expectation is indeed borne out 
in practical terms: Even though strictly speaking there 
cannot be a first-order transition with any amount of 
disorder, ESI the smearing of the transition is very small 
and the observable effect is just a suppression of the tri¬ 
critical temperature, with the quantum phase transition 
remaining first orderHowever, for a threshold value of 
the disorder strength the tricritical temperature reaches 
zero, the quantum phase transition becomes second or¬ 
der, and above this threshold it remains second order, 
albeit with an unusual critical behavior. The evolution 
of the phase diagram is shown in Fig. [3] A renormalized 
mean-field theory leads to unusual critical exponents, and 
a renormalization-group analysis of the fluctuations re¬ 
veals that there are marginal operators in all dimensions 


Figure 3: Evolution of the phase diagram in temperature- 
control parameter-field ( T-r-h ) space with increasing disor¬ 
der. For weak disorder the phase diagram is a schematic 
version of the experimentally observed one shown in Fig. [2] 
Increasing disorder suppresses the tricritical point (TCP). At 
a critical disorder the TCP merges with the quantum wing- 
critical points (QWCPs), and a quantum critical point (QCP) 
in zero field emerges. For disorder strengths close to the 
threshold the asymptotic quantum critical behavior is con¬ 
fined to a very small region (red circle), and the observable 
critical behavior (blue circle) is controlled by Hertz’s fixed 
point. With further increasing disorder the instability of 
Hertz’s fixed point becomes apparent farther away from the 
transition, and the region controlled by the true critical fixed 
point grows. 


2 < d < 4, which leads to log-normal terms that multiply 
the usual critical power laws.E^HHl This modification of 
scaling mimics, in a sizable pre-asymptotic region, power 
laws with effective e xpo nents that are quite different from 
the asymptotic ones PS Furthermore, the behavior associ¬ 
ated with Hertz’s unstable fixed point can be observable 
in substantial regimes in parameter space, as we will show 
below. We therefore list in Table |T] and discuss below, 
the critical behavior at Hertz’s fixed point as well as that 
associated with the physical critical fixed point and its 
pre-asymptotic region. 


1. Hertz’s fixed point 

As mentioned above, Hertz’s fixed point in the disor¬ 
dered case is never a physical fixed point; it is unsta¬ 
ble with respect to the same physical processes as in the 
clean case, viz., coupling of the magnetization to soft 
fermionic excitations. For disordered systems, the lat¬ 
ter are diffusive, and their effects are small in the limit 
kpt 1, with kp the Fermi wave number and t the elas¬ 
tic mean-free path. In addition, it is destabilized by a 
random-mass term in the LGW theory for the magnetic 
order parameter t Allifl This effect is small as long as the 
fluctuations of the distance from criticality, Sr , are small 
compared to r itself: Sr oc (A/^) 1 / 2 oc (Ar^) 1 ^ 2 <C r. 
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Here A is the disorder strength, and we assume that the 
fluctuation is proportional to the inverse square root of 
a correlation volume £ d , with £ the correlation length. 
With mean-field exponents in d = 3, this condition be¬ 
comes A 2 <C r. With A = 1 /kp£ this is a very weak 
constraint for small disorder, kpi 1, and Hertz’s mean- 
field description remains valid in a large parameter range. 
In d = 3, the condition becomes £ -C t. In Fig. [3] the 
region controlled by Hertz’s fixed point is schematically 
indicated by the blue region around the quantum criti¬ 
cal point, and the asymptotic critical region by the red 
circle inside the blue region. We conclude that in many 
systems one expects sizable regions in the phase diagram 
where Hertz’s fixed point yields the observable behavior, 
and only very close to the quantum critical point does 
the behavior cross over to the asymptotic one. This ex¬ 
perimental relevance of Hertz’s fixed point in disordered 
systems has not been appreciated before, and we there¬ 
fore discuss it here. 

Table |T] lists the critical exponents associated with 
Hertz’s fixed point in the presence of disorder. The static 
exponents have their mean-field values for all d > 0, and 
the dynamical exponent z c = 4 results from the diffusive 
dynamics of the conduction electrons, see above. The 
temperature dependence of the observables, expressed 
by the T-exponents, is obtained by repeating Millis’s 
analysis with obvious modifications; the most important 
one being that the scale dimension of the DIV is now 
[u] = —d, rather than [u] = — (d — 1) in the clean case. 

The system is above its upper critical dimensionality 
for all d > 0, and the exponent relations that depend 
on strong scaling therefore break down, while those that 
rely only on weak scaling hold, in perfect analogy to the 
clean case in d > 1. The discussion of the exponents 
a and o.t in Sec. |III A 1| also carries over with obvious 
modifications. At the upper critical dimension d+ = 0 
hyperscaling holds as expected, and there is only one 
dynamical exponent, z m = z c = 4, as is the case in clean 
systems in d = 1. An important exponent relation that 
has no clean analog is the rigorous inequality v > 2/ dW^ 
This is violated for all d < 4, which by itself implies that 
the fixed point cannot be stable. 

We next consider the temperature dependence of the 
conductivity. The arguments that led to Eq. (3.2) are 


easily modified to apply to the disordered case. Since 
the scale dimension of the DIV is now [it] = —d, the 
effective scale dimension of Act is [Act] = — (d + 4). In 
addition, the dynamical exponent is now z c = 4, which 
yields 


Aer(r, T) = b d+4 F a (r ft 1 /", T b 4 ) . 


(3.6) 


This immediately yields St = — {d + 4)/4 as shown in 
Table |T] Here we have assumed that the backscattering 
factor is still present, which requires that the residual re¬ 
sistivity is small compared to the temperature-dependent 
part; see the discussion below. The exponent s does not 
exist in this regime for the same reason as in the clean 


The disorder dependence of the temperature- 
dependent conductivity can also be obtained from 
scaling. Let us start from Eq. (3.21 with z c = 3. 


Dropping the meaningless r-dependence, and adding the 
dependence of Act on the elastic mean-free path £, we 
have 


A a{T,£) = b d+2 F a {Tb 3 ,£b~ 1 ) 


(3.7) 


The scaling function F a now must have the following 
prop ertie s: (1) F a (x,y —>• oo) = F^ ean (x), so we recover 
Eq. (3.2). (2) For t/b —» 0 the prefactor must change 
from b a+2 to b d+4 , and the temperature scale f acto r must 
change from b 3 to b 4 in order to recover Eq. (3.61. This 
is achieved by F a (x,y —>■ 0) = F dlsordered (x /y) / y 2 . We 
thus find the following generalization of Eq. (3.6): 


Act(T) = b d+4 £~ 2 F a ((T/£) b 4 ) . (3.8) 


This yields 


Act(T) a £~ (4 " d)/4 T- (d+4)/4 , (3.9) 


which yields the disorder dependence of the prefactor in 
addition to the exponent st ■ We have checked this result 
by means of an explicit calculation based on the Kubo 
formula for the conductivity, and have found agreement. 

The above results for the conductivity depend on var¬ 
ious assumptions that limit their regime of validity and 
require a discussion. Firstly, we have assumed diffu¬ 
sive electron dynamics, which requires TV < 1, with 
r = £/vp the elastic mean-free time. Secondly, we have 
assumed that the contribution from Eq. (3.91, A p oc 
£(4—d)/4 y(d+4)/4, w } 1 i c ] 1 j s contribution from small 
wave numbers in any transport theory, is larger than the 
clean contribution, A p oc T^ d+2 ^ 3 . Thirdly, we have as¬ 
sumed that the backscattering factor is still present. This 
requires that the residual resistivity po be small compared 
to the temperature-dependent contribution p(T). The 
last two assumptions both imply that the temperature 
must be large compared to a disorder-dependent energy 
scale, and the result (3.91 will thus be valid only in a tern- 


case. 


perature window. In a simple model for a metallic ferro- 
magnet, where the Fermi energy is the only microscopic 
energy scale, this window does not exist, since the last 
requirement leads to an unrealistically large lower bound 
for the temperature. This is misleading, however. In any 
real ferromagnetic material a complicated band struc¬ 
ture leads to multiple microscopic energy scales, some 
of which are strongly renormalized downward from the 
Fermi energy of a nearly-free-electron model. This is es¬ 
pecially true for the low-Curie-temperature ferromagnets 
that are good candidates for observing a quantum ferro¬ 
magnetic transition. Different factors of the temperature 
in the above arguments will be normalized by different 
microscopic scales, and generically one expects the rele¬ 
vant temperature window to exist. This is certainly true 
empirically, as the temperature-dependent resistivity is 
observed to dominate the residual resistivity for all but 
the lowest temperatures in many materials] 47 ! 
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In the regime discussed above the exponent s as de¬ 
fined in Eq. (A 6 l does not exist for the same reason as 
in the clean case, and the low-temperature behavior of 
the conductivity away from criticality is again given by 
Eq. (3.31. Remarkably, the r-dependence of the prefactor 
of the T~ 2 dependence is the same in the clean and dirty 
cases. 


2. The physical fixed point 


Asymptotically close to the quantum critical point the 
behavior is described by a Gaussian fixed p oint w ith 
marginal operators in all dimensions 2 < d < 4 j 4 :j l 44 l an d 
with increasing disorder the region in parameter space 
that is controlled by this fixed point grows. The marginal 
operators result in critical behavior that is not given by 
pure power laws. For example, the magnetization at 
T = 0 in d = 3 asymptotically behaves as 

m(r,T = 0) a (—r ) 2 [g(ln(—1/r ))] 2 . (3.10a) 

Here the function g is asymptotically log-normal, 

g(x —> oo) oc e cx , (3.10b) 

with c a constant; see Ref. 3H for details. In a large 
pre-asymptotic region these multiplicative corrections to 
scaling mimic power laws that span several decades in r 
or T.ESFor instance, the magnetization obeys an effective 
homogeneity law 

m(r,T,h) = b~ 2 F m {rb d - 2+x ,Tb 2+x ,hb d+x ) . (3.11) 


A is an effective exponent that depends on d and r and 
goes to zero for r —> 0 ; however, it is approximately 
constant over a large r-range. For d = 3, A « 2/3 for 
0-001 l r l 0.1. As an example, we show in Fig. [4] 
the divergence of the specific-heat coefficient, l/q(r, T = 
0) oc |r|“. The asymptotic value A = 0 needs to be inter¬ 
preted as signaling the presence of the log-normal terms. 
The resulting power laws are listed in the two middle 
columns in Table [I] Since the system is below an upper 
critical dimension for 2 < d < 4, all exponent relations 
hold, including the hyperscaling relations. This holds 
even for the pre-asymptotic effective exponents, with the 
exception of the generalized Harris criterion, Eq. (Cl I, 
which holds only asymptotically. 


IV. SUMMARY, AND DISCUSSION 
A. Summary 

In summary, we have presented two conceptual ad¬ 
vances for quantum phase transitions: First, we have an¬ 
alyzed relations between critical exponents, also known 
as scaling relations. Compared to classical transitions, 
there are additional exponents to consider that describe 


TT 

r 
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Figure 4: The inverse of the specific-heat coefficient 7 nor¬ 
malized by the Fermi-liquid value 7 0 as a function of r in the 
ordered phase for d = 3. The inset demonstrates the effec¬ 
tive power law 7 oc |r| _ “ with a = A/(l + A) ~ 0.4. Figure 
adapted from Ref. 1451 


the scaling behavior of observables with respect to tem¬ 
perature, which is independent from the scaling behav¬ 
ior with respect to the control parameter. The expo¬ 
nents describing the latter (“r-exponents”) are related to 
the ones describing the former (“T-exponents”) by means 
of a dynamical critical exponent; however, care must be 
taken since at many quantum phase transitions there is 
more than one dynamical exponent. We have shown that 
the Widom and Fisher equalities hold for both the r- 
exponents and the T-exponents under weak assumptions 
that are fulfilled at most quantum (and classical) phase 
transitions. The same is true for the Essam-Fisher equal¬ 
ity relating the r-exponents. Additional relations, which 
are akin to the classical hyperscaling relation between the 
correlation-length exponent v and the specific-heat expo¬ 
nent a, require stronger assumptions that break down if 
the system is above an upper critical dimension. We 
have also generalized the rigorous Rushbrooke inequality 
to the quantum case, and we have shown that at quan¬ 
tum phase transitions it is crucial to distinguish between 
exponents describing the critical behavior of the specific- 
heat coefficient, which we denote by a and cxt, and the 
exponents a and ar that govern the control-parameter 
susceptibility; at a classical transition, a and a are the 
same. We have also discussed Wegner’s equality that re¬ 
lates critical exponents for the electrical conductivity and 
the correlation-length exponent v, and have discussed the 
conditions under which it is valid. 

Second, we have generalized the concepts of Fisher and 
Berker related to scaling at classical first-order transi¬ 
tions to quantum phase transitions. The scaling concepts 
all carry over, with the temperature playing the role of 
additional dimensions. Again, the presence of multiple 
time scales leads to complications that need to be dealt 
with carefully. With the proper choice of dynamical ex¬ 
ponents all scaling relations, including the hyperscaling 
ones, hold under very weak assumptions. This reflects 
the fact that at a first-order transition there are no dan- 
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gerous irrelevant variables. 

We then have applied these concepts to the case of the 
ferromagnetic quantum phase transition as a specific ex¬ 
ample. This transition is well suited for this purpose, 
since it shows very different behavior in the presence or 
absence of quenched disorder, respectively, and in either 
case there are interesting and experimentally relevant 
crossover phenomena that are governed by different fixed 
points, some of which are above and some of which are 
at their upper critical dimension in the physically most 
interesting dimensions d = 2,3. 


B. Discussion 

We conclude with some discussion points that augment 
remarks already made in the body of the paper. 


1. General points 

Number of independent exponents: There are seven 
r-exponents for thermodynamic quantities: a, a , /3, 7 , 
<5, 77 , and v. These are related by three independent 
weak-scaling expone nt rela tions (Widom, Eq. ( 2.13a l, 
Essam-Fisher, Eq. (2.15aI, and Fisher, Eq. (2.20al), 
and two independent strong-scaling ones (a-hyperscaling, 
Eq. (2.23aI, and a hyperscaling, Eq. (2.24a). In the pres¬ 


ence of strong scaling one thus has two independent static 
r-exponents. In addition, there are five T-exponents: ar, 
ax, Pt, and 7 t that are constrained by two indepen¬ 
dent weak-scaling exponent relations (Eqs. ( |2.13b[[Z20bl) 


plus three independent strong-scaling ones (Eqs. (2.15b 


2.23b 2.24b l. In the presence of strong scaling, we thus 


have two independent static exponents, plus the dynamic 
ones. In the absence of strong scaling, there are ei¬ 
ther four or five independent static exponents plus the 
dynamic ones, depending on whether the Essam-Fisher 
equality for the T-exponents holds. 

Observability of the exponent Pt The order parame¬ 
ter is nonzero only in the ordered phase, i.e., for — r > 
rpi/vzm (ignoring constant proportionality factors). For 
-r > T x l VZm one observes static scaling with small tem¬ 
perature corrections, and for — r < T 1 / VZm the scaling 
function vanishes identically. The exponent (3t therefore 
cannot be observed via the T-dependence of m at r = 0. 
Instead, consider the general weak-scaling homogeneity 
law for m 

m(r ,T) = b~ f) / v § m {rb 1/v , Tb f>,v ^ T ) , (4.1a) 


which follows, e.g., by differentiating Eq. (2.6 1 with re¬ 
spect to h and using Eqs. (2.13a) and (2.15a). This can 
be written 


m(r,T) = T^$ m (r/T^,l) 
= r^ m (l,T/r^) , 


(4.1b) 


which defines the exponent Pt- Note that 4> m (a:, 1) has a 
zero for some value x = x c < 0 that determines the phase 
boundary, and that & m {x > x c , 1) = 0. For the temper¬ 
ature derivative of the order parameter this implies 


dm 

dT 


(r,t) = T^- 1 l> m (r/T^,l) , 


(4.2) 


where $ m ( x,y ) = d y 4> m (a:, y ). This has been used to 
derive Eq. ( 2.22a| ) from ( |2.21 ). 

As an illustration, consider the third panel of Fig. [T] 
again. The data imply m 2 oc const. — T 4//3 , where the 
constant depends on r, but the prefactor of the T 4//3 does 


not. When interpreted by means of Eq. (4.1b I this im¬ 


plies p = 1/2 and Pt/P = 4/3. The latter ratio is equal 
to 1 /vz m , which is consistent with the observed shape of 
the phase boundary in the first panel of the figure. All of 
this is consistent with the discussion in Sec. IIII A II An¬ 
other way to measure Pt is to observe the T-dependence 
of the susceptibility at criticality, which yields 7 ^, see the 
second panel in the figure, and use Eq. (2.13b) in conjunc¬ 


tion with an independent determination of the exponent 
S. We finally mention that for asymptotically small m at 
any nonzero temperature one crosses over to the classical 
critical behavior characterized by the classical value of 
the exponent p. 

Coupling of statics and dynamics: Even though the 
statics and the dynamics are coupled in quantum statis¬ 
tical mechanics, the Widom, Essam-Fisher, and Fisher 
exponent equalities do not involve the dynamical expo¬ 
nents. Even at T = 0 there thus are exponents relations 
that do not reflect this coupling. 

Relevance of quenched disorder: Since the generalized 
Harris criterion by Chayes et al., Eq. (Cl I, is rigorous 


it provides a necessary condition for the stability of a 
critical fixed point in the presence of quenched disorder. 
For instance, it immediately tells us that Hertz’s fixed 
point in disordered systems with its mean-field value v = 
1/2 cannot be stable for d < 4, although it does not 
provide any hints as to what the fixed point is unstable 
against, or what replaces it asymptotically. Similarly, 
it shows that a first-order transition is strictly speaking 
impossible in a disordered system, since it would result 
in a correlation-length exponent v = \/{d + z m ) < 2 /d. 

Significance of a = 1 at a first-order QPT: A first- 
order classical phase transition is characterized by the 
appearance of a latent heat, i.e., a discontinuity of the en¬ 
tropy as a function of the temperature, and a (5-function 
contribution to the specific heat. The corresponding 
physical phenomenon at a first-order QPT is a discon¬ 
tinuity of the derivative of the free energy with respect 
to the control parameter, and a d-function contribution 
to the second derivative. As an example, consider a QPT 
where the control parameter is the hydrostatic pressure p. 
With g the Gibbs free energy density, the compressibility 
k = —d 2 g/dp 2 will thus have a (5-function contribution. 
That is, at any pressure-driven QPT the system will dis¬ 
play a “latent-volume” effect, i.e., the system volume will 
change spontaneously and discontinuously at the critical 
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pressure. 

Choice of the control parameter: While deriving a 
LGW theory, the easiest choice for the control parameter 
or mass term r may be a linear combination of some non- 
thermal parameter and the temperature, or some power 
of the temperature. This is not a good choice for r for 
at least two reasons: (1) In order to measure, e.g., the 
exponent /3 it is necessary to keep T = 0 in addition to 
keeping the conjugate field equal to zero. Moving on a 
path that is not in the T = 0 plane is akin to deviating 
from the critical isochore at a classical liquid-gas tran¬ 
sition. (2) The such-obtained temperature dependence 
of the phase diagram may not be the leading one due 
to DIVs. This is what happens, for instance, in Hertz 
theory. 

Applicability of scaling theory: We have restricted 
ourselves to phase transitions for which a local order pa¬ 
rameter exists. We note, however, that the existence or 
otherwise of a local order parameter can be a matter of 
how the theory is formulated. In the case of the fer¬ 
romagnetic quantum phase transition in metals that we 
have used as an example in Sec. |III| it is crucial that one 
treats the fermionic soft modes that couple to the order 
parameter explicitly, as in Ref. 1431 if one integrates them 
out to formulate a theory entirely in terms of the order 
parameter a local description is not possible.® We also 
note that some theories of “exotic” quantum phase tran¬ 
sitions that cannot be cast into the language of a simple 
Landau theory are structurally very similar to the ferro¬ 
magnetic quantum phase transition problem in that they 
couple an order-parameter field to a gauge field,® or to 

fermions.® 

More generally, we note that the scaling arguments we 
have employed are extremely general and hinge only on 
the existence of a phase transition with power-law critical 
behavior. Given this, the free energy will obey a general¬ 
ized homogeneity law irrespective of whether or not the 
transition allows for a traditional Landau description or 
is of a more exotic nature. 


2. Points related to Hertz’s fixed point 

Multiple time scales in Hertz theory: We come back to 
the identification of the dynamical critical exponent z m 
at the beginning of Sec. |III A 1| It is important to dis¬ 
tinguish between theories that truly have multiple time 
scales that belong to different types of excitations, and 
theories where a single time scale may or may not get 
modified by the effects of DIVs, depending on the context, 
which effectively leads to multiple time scales. Hertz’s 
theory, both for the clean and the disordered case, be¬ 
longs to the latter category. The time scale with z m = 1 
in the clean case, or z m = 2 in the disordered one, at the 
physical fixed point (see Table [I]) is missing in Hertz the¬ 
ory, and the effective z m listed for Hertz’s fixed point in 
Table |T] is due to a DIV modifying the sole critical z = 3 
(clean) or z = 4 (disordered). Note that if a 0 = 1 were 


present in the clean theory, it would be smaller than the 
z m = 6 /(d + 1) in Hertz-Millis theory, and thus would 
play the role of z m according the arguments in the con¬ 
text of Eqs. ( |2.8| |2.11 1 . This is another indication that 
Hertz’s fixed point cannot be stable. 

Numerical values of exponents: The values of the 
commonly measure exponents dp, 7 t, and ftp are rather 
similar at the clean and dirty Hertz fixed points, respec¬ 
tively, in d = 3. The same is true for the scaling of the 
Curie temperature with the control parameter, which is 
Tc oc r 3 / 4 in the clean case, and Tq oc r 4//5 in the dis¬ 
ordered one. This needs to be kept in mind when inter¬ 
preting experiments. 

Also of interest is the value of the conductivity expo¬ 
nent st = —7/4 at the disordered Hertz fixed point in 
d = 3. This result is a good candidate for explaining the 
fact that the resistivity is commonly observed to vary 
as T x with x > 1 near a ferromagnetic quantum critical 
point, as the crossover to the asymptotic critical behavior 
is expected to occur only at very low temperatures. 


3. Points related to the physical quantum ferromagnetic 
fixed point 


Significance of z c = d: In Sec. Ill A 2 we showed that 
z c = d follows from the first-order nature of the quan¬ 
tum phase transition. Combined with the origin of the 
dynamical exponent in a LGW theory, where z c arises 
from a combination of the Landau-damping term with 
the spin susceptibility in the paramagnetic phase, this 
implies that in a clean Fermi liquid the leading non- 
analytic wave-number dependence must scale as k d ~ 4 < 
The connection between these seemingly unrelated re¬ 
sults is scaling, see Sec. |III A2| and Ref. [SD] 

In this context we also note that at any first-order 
quantum phase transition in a fermionic system one ex¬ 
pects the specific-heat coefficient to be discontinuous. 
This implies a = dtp = 0, which in turn requires z c = d, 
see Eq. (2.33g). 


Scale dimension of the conjugate field: We saw in 
Sec. |II C| that a dimensionless order parameter implies 
a scale dimension [h] = d + z m for the conjugate field. 
This illustrates the fact that in general it is not possible 
to simply relate the scaling of h to the scaling of the en¬ 
ergy or temperature. For instance, one might argue that 
the scale dimension of h should be equal to z c , since h 
determines the Zeeman energy. While this is sometimes 
true (for instance, at the disordered critical fixed point), 
it is not true in general. In this context we also note that 
the Rushbrooke inequalities, Eqs. (2.221, do not depend 


on [h], as [h] cancels between the contributions 2fix and 
7 t (or 2/3 and 7 ), respectively, on the left-hand side. 

Logarithmic corrections to scaling in a range of dimen¬ 
sions: The theory of Refs. l l.'lllTI yields logarithmic cor¬ 
rections to scaling (in the sense of log-normal terms mul¬ 
tiplying power laws) in an entire range of dimensions, 
2 < d < 4. This is unusual; more commonly logarith- 
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mic corrections to scaling occur only in a specific dimen¬ 
sion. This can be understood within Wegner’d^ classifi¬ 
cation of logarithmic corrections, as has been discussed 
in Ref. [3Tl The dynamical critical exponent z c = d leads 
to operators that are marginal in a range of dimensions, 
and this in turn causes the logarithmic terms. We also 
note that the mathematical problem of determining the 
critical behavior for the theory of Refs. 1-13 ITT! was solved 
exactly in Ref. [221 although the physical interpretation 
was unclear at that time. See also Ref. f53l 

Scaling of the Griineisen parameter The Griineisen 
parameter T = (— ds/dp)/T(ds/dT ), with s the entropy 
density, is defined as the ratio of the thermal expansion 
coefficient (dV/dT) p ^/V = —{ds/dp)T,N and the spe¬ 
cific heat c = Tds/dT. It was shown in Ref. [5U that 
at a pressure-tuned QPT with a single dynamical expo¬ 
nent z the Griineisen parameter diverges as T oc T -1 / 1 ' 2 ; 
this is readily confirmed by using Eq. (2.51. With two 
dynamical exponents Z\ > Z 2 we find from Eq. (2.8) 


T oc In particular, at the physical fixed point in 

disordered metallic ferromagnets we have T oc T -1 / 3 in 
d = 3. Another useful observation is that, according to 
its definition, T scales as T ~ 1/p ~ 1 /r at a pressure- 
tuned QPT. Since r ~ T 1 /" 2 at a critical point, this is 
equivalent to the result given above. At the first-order 
transition in clean systems, T ~ 1 /r reflects the same 
(5-function contribution that was discussed for the com¬ 
pressibility in Sec. |IVB 1| This is also apparent from 
the fact that s ~ 7 T, and hence T ~ d'y/dp, and the 
specific-heat coefficient 7 has a discontinuity across the 
first-order transition. 

The behavior at the clean Hertz fixed point has been 
discussed in Ref. [54] Up to logarithmic corrections, T ~ 
r _1 ~ T -2 / 3 in agreement with the above arguments and 
v = 1/2, z c = 3 from Table [I] At the disordered Hertz 
fixed point the corresponding result is T ~ r _1 ~ T -1 / 2 . 
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Appendix A: Definitions of critical exponents 

Let T be the temperature, h the field conjugate to 
the order parameter, r the control parameter, i.e., the 
dimensionless distance from criticality at T = h = 0, and 
/ a suitable free-energy density. Consider the correlation 
length £, the order parameter m = df/dh, the order- 
parameter susceptibility Xm = d 2 f/dh 2 , and the specific- 
heat coefficient 7 = c/T = —d 2 f/dT 2 as functions of r, 
T, and h, and the susceptibility also as a function of the 
wave number k. Further consider the susceptibility \r = 
— d 2 f/dr 2 , which we refer to as the control-parameter 
susceptibility. (More precisely, \r is the susceptibility 
of the thermodynamic quantity whose conjugate field is 


the control parameter.) We define critical exponents at 
a quantum phase transition as follows. 

Correlation length: 

£(r -A 0,T = 0) oc \r\~ v 
Order parameter: 


£(r = 0, T -a 0) ocT~ 


(Al) 


m(r — t 0, T = 0, h = 0) oc (— r)^ 
m(r = 0, T = 0, h —► 0) oc h 1 ^ 5 , 
m(r = 0,T-a 0,h = 0) oc T /3t . 


(A2) 


The last definition is purely formal: Since the magneti¬ 
zation is nonzero only in the ordered phase, the T^ T in 
the last line has a zero prefactor. See the discussion in 


Sec. IV B 1 of how to interpret the exponent /?t- 
Order-parameter susceptibility: 

Xm( r —t 0, T = 0; fc = 0) a |rj -7 , 

Xm(r = 0, T -A 0; fc = 0) cx T ~, 

X m{r = 0,T = 0,fc—>0) cx 1/k 2 ~ v . 

Specific-heat coefficient: 

7 (r —>• 0, T = 0) oc |r| _ “ , 7 (r = 0, T —> 0) oc T~ 
Control-parameter susceptiblity 


(A3) 


(A4) 


Xr(r-+0,T = 0)oc |r| _Q , Xr (r = 0, T -a 0) oc T~ aT . 

(A5) 

v , /?, 7 , 8, and 77 are defined in analogy to the corre¬ 
sponding exponents at a classical phase transition . 1 In 
the main text we refer to these exponents, and also to a 
and a, as the /’-exponents. The definition of a deviates 
from the one of the classical exponent customarily de¬ 
noted by a, which is defined in terms of the specific heat 
rather than the specific-heat coefficient. This is necessary 
in order to factor out the factor of T in the relation be¬ 
tween the specific heat and the specific-heat coefficient,^ 
which makes no difference at a thermal phase transition, 
but goes to zero at a QCP. For instance, the thermody¬ 
namic identity that underlies the Rushbrooke inequality, 
Eq. (B5l, has no explicit T-dependence only if it is for¬ 


mulated in terms of specific-heat coefficients rather than 
the specific heats. At a classical phase transition, a coin¬ 
cides with a. Our definition of a at a QPT is analogous 
to the definition of the classical exponent a from a scal¬ 
ing point of view. However, the physical interpretation of 
a and the susceptibility \r depends on the nature of the 
control parameter. For instance, if the control parame¬ 
ter is hydrostatic pressure, and / is the Gibbs free energy 
density, then \r is proportional to the compressibility of 
the system. At a thermal transition, where r oc T — T C1 
Xr is identical with the specific-heat coefficient 7 , and a 
has its usual meaning, or, dr> vt> Pt, and 7 t reflect the 
fact that a QPT can be approached either in the T = 0 
plane, or from T > 0. These exponents are referred to as 
the T-exponents in the main text. 










18 


We finally define exponents s and st that describe the 
behavior of the electrical conductivity a at the QPT. 
With Act the scaling part of a, these are defined as fol¬ 
lows: 

Electrical conductivity: 

Acr(r-S>0,T = 0)oc \r\ s , Aa(r = 0, T -A 0) oc T ST . 

(A 6 ) 

The exponent s obviously makes sense only in systems 
with quenched disorder. 


Appendix B: Classical scaling relations 

For the convenience of the reader we recall some well- 
known classical exponent relations (see, e.g., Refs, [mil : 


t-H 

1 

'-o 

<=CL 

II 

(Widom) 

(Bla) 

a + 2/3 + 7 = 2 

(Essam-Fisher) 

(Bib) 

II 

to 

1 

(Fisher) 

(Blc) 

v d = 2 — a 

(“hyperscaling”) 

(Bid) 


where a is the ordinary specific-heat exponent. The first 
two equalities follow from a weak scaling assumption, in 
the sense of Sec. |II A| for the singular part / of the free- 
energy density. This has the effects of dangerous irrel¬ 
evant variables, if any, built in and can be written, for 
instance, as 

f(r,h) = b-^- a ^ v ^ f [rb xlv ,hb^ lv ) , (B2) 


with $/ a scaling function. The Fisher scaling rela¬ 
tion, Eq. ( Blc I, requires an additional (weak) scaling 
assumption,® namely, that scaling works for correlation 
functions as well as for thermodynamic quantities. If we 
use a homogeneity law for the wave-number dependent 


order-parameter susceptibility, 

x(r,k) = b^ 1 2 3 4 5 6 '$ x (rb 1 ' I ',kb) , (B3) 

and put r = 0 we obtain Eq. ( |Blc| ). The hyperscaling re¬ 
lation requires yet another assumption,® which is a strong 
scaling assumption in the sense of Sec. |II A| and which is 
equivalent to saying that the scale dimension of the free 
energy density is equal to d: 

f(r,h) = b- d ^ f {rb 1/v ,hb 0S/v ) . (B4) 


If hyperscaling holds, the six exponents a, /3, 7 , <5, rj, and 
v are constrained by the four relations in Eqs. (B1), and 
only two of the exponents are independent. 

We also list the Rushbrooke inequality 


a + 2/3 + 7 > 2 , (B5) 


which holds rigorously, as it only depends on thermody¬ 
namic stability conditions. There are other inequalities 
that depend on various assumptions ! 1 

Appendix C: The Harris criterion 


Harris considered the effects of quenched disorder on 
classical critical points. He found that the fixed point de¬ 
scribing the transition in the clean system is stable with 
respect to disorder, i.e., the critical behavior is unaffected 
by the disorder, as long as the specific-heat exponent is 
negative, a < 0.^ This statement is referred to as the 
Harris criterion. Chayes et aP^3 generalized Harris’s ar¬ 
gument and made it rigorous to show that, under very 
general conditions, the correlation-length exponent at a 
critical point in a disordered system, classical or quan¬ 
tum, must obey 

v > 2/d (Chayes-Chayes-Fisher-Spencer) . (Cl) 
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